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Abstract Many of the challenges of scaling quantum computer hardware lie at the 
interface between the qubits and the classical control signals used to manipulate them. 
Modular ion trap quantum computer architectures address scalability by constructing 
individual quantum processors interconnected via a network of quantum communi¬ 
cation channels. Successful operation of such quantum hardware requires a fully pro¬ 
grammable classical control system capable of frequency stabilizing the continuous 
wave lasers necessary for trapping and cooling the ion qubits, stabilizing the optical 
frequency combs used to drive logic gate operations on the ion qubits, providing a 
large number of analog voltage sources to drive the trap electrodes, and a scheme for 
maintaining phase coherence among all the controllers that manipulate the qubits. In 
this work, we describe scalable solutions to these hardware development challenges. 
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1 Introduction 

Trapped ions exhibit many favorable properties for use as qubits |[Il|2l, such as long 
coherence times m, near perfect state preparation and readout mi, and high fi- 
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delity qubit gate operations im). The Modular Universal Scalable Ion-trap Quan¬ 
tum Computer (MUSIQC) platform provides a scalable architecture where registers 
of ion qubits can be trapped in scalable surface electrode ion traps ll9lfT0l . Each quan¬ 
tum register must be capable of performing arbitrary quantum logic operations among 
the qubits in the register, as well as supporting the optical interface needed to com¬ 
municate qubits between quantum registers through a photonic network CHEl. 

Any ion trap quantum computer, and specifically the MUSIQC architecture ifTTl 
El, requires a wide range of classical hardware to drive the quantum operations 
needed for the computation. Scaling up from a small laboratory experiment to a large- 
scale processor will require scalable classical resources in addition to qubits. Here we 
present scalable hardware to perform some of the classical tasks required by an ion 
trap quantum computer. These tasks include: continuous wave (CW) laser frequency 
locking, frequency stabilization of pulsed laser beams for quantum logic gates, and 
low-noise digital-to-analog converter (DAC) systems for trapping and shuttling ions. 
We designed each of these systems to be inexpensive, compact, and extensible such 
that they can grow or be replicated to accommodate a much larger experiment. The 
techniques are developed for ^^^Yb^ hyperfine qubit levels separated by approxi¬ 
mately 12.6 GHz, where resonant laser beams at 369.5 nm and off-resonant Raman 
beams detuned 14-33 THz from the 1/2 to ^Pi /2 levels are used to manipulate the 
qubits. Many of these techniques can be extended to hyperfine qubits in other ion 
species. 


2 Offset frequency lock for CW laser stabilization 

The narrow spectral linewidths (< lOMHz) that are typical of atomic qubits and neu¬ 
tral atom experiments require lasers with high frequency specificity and stability. 
Current state-of-the-art laser systems that rely on the control of operating current and 
temperature are unable to sufficiently stabilize the frequencies of commercial laser 
diodes to meet these needs, so active locking to another frequency reference is of¬ 
ten used. Commercially available vapor cells allow for easy stabilization of lasers 
addressing neutral atom transitions, but such a simple and inexpensive solution is 
not available to lasers addressing ion transitions (by virtue of their ionic nature), and 
other methods must be employed. 

One approach to frequency locking an external cavity diode laser (ECDL) is to 
use a second laser, locked to a vapor cell at frequency /vapor (e.g. 780 nm for Rb 
vapor), and an optical cavity. The optical cavity is designed to support both /vapor 
and light at the desired ECDL frequency (e.g. 369.5 nm for ^^^Yb^ ions). The op¬ 
tical cavity length is stabilized to /vapor while the cavity error signal is used to pro¬ 
vide feedback to the piezo-controlled grating inside the ECDL, effectively transfer¬ 
ring the vapor cell frequency reference to the ECDL. A second approach, where the 
ECDL is stabilized to the transitions of an ion plasma (for example, in a hollow- 
cathode lamp) can also provide the frequency stability necessary for trapped ion ex¬ 
periments ca. While both approaches can adequately stabilize the frequency of an 
ECDL for trapped ion experiments, scaling is a concern as they are expensive and 
cumbersome to implement, especially for a complicated experiment involving multi- 
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Fig. 1: Schematic of the offset lock setup. A portion of the output power of an ab¬ 
solutely frequency stabilized master ECDL (M) is combined with a number of slave 
ECDLs (S) on a photodiode (PD). The measured beat frequency is then divided by a 
1 /N prescaler {N = 10, FPS-10-12, RF-Bay) and compared to the setpoint frequency 
output of a direct digital synthesizer (DDS) (/dds) by a frequency/phase detector. 
The resulting error signal is filtered by a proportional-integral (PI) controller and fed 
back to the piezo-cotrolled grating of the slave ECDL to equalize the beat and set- 
point frequencies making the frequency offset between the master and slave laser 
^ X /dds- 


pie lasers. Here we present a scheme that transfers the lock of a single master laser to 
multiple slave lasers with minimal hardware, complexity, and cost. 

In this scheme (shown in Fig. [T]), a portion of the slave’s output and a small 
amount of light from the frequency stabilized master laser are matched in mode shape 
and polarization, and overlapped on a photodiode (PD). The overlapped electric fields 
result in a beat note at the difference frequency between the master and slave lasers, 
which is detected by PD. The range of detectable beat note frequencies is limited 
by the PD bandwidth and is typically a few tens of MHz to several GHz. The beat 
note is amplified and compared to a reference signal (generated by a direct digital 
synthesizer (DDS) in our case) using a phase/frequency detection circuit. After ap¬ 
propriate filtering, the output of the circuit can be used as an error signal to provide 
feedback to the piezo-tunable grating and/or current of the slave ECDL using stan¬ 
dard Pound-Drever-Hall control techniques l(T4l . The frequency difference between 
the master and slave lasers is digitally programmable by the DDS frequency to be 
N X /dds using a 1/A prescaler, as shown in Fig.[^ 

This solution is highly scalable as many slave lasers can be locked to a single 
master laser. The number of slave lasers that can be locked to a single master laser is 
limited by the total power of the master laser. In our system, a 369.5 nm ECDL, with 
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Fig. 2: Allan deviation of the beat note frequency as a function of averaging window 
time for a stabilized 369.5 nm ECDL. The Allan deviation shows 1/v^ scaling and 
is indicative of white frequency noise. Inset: The difference between the locked beat 
note frequency and target frequency as a function of time for a 1 second averaging 
window as measured by a frequency counter. Data shown is one data point per second 
for one hour. 


an output power of 3.4 mW, serves as the master laser which is frequency stabilized 
to a hollow-cathode lamp ca. From the master laser, 600 /iW of power is used to 
lock 3 slave lasers at different offset frequencies (determined by the DDS frequency). 
We are currently capable of locking up to 8 slave lasers to a single master laser, 
which is more than enough to meet our current experimental demands. If necessary, 
this number could be pushed higher with further improvements in mode-matching, or 
with more specialized equipment. The beat note is detected by a high-speed response 
(2 GHz) photodiode (S9055, Hamamatsu), thus we are able to stabilize the master 
laser to the ^^^Yb^ and have slave lasers stabilized to either the Yb^ and ^^'^Yb^ 
resonances as necessary. From the Allan deviation of the beat note frequency for a 
slave 369.5 nm ECDL, we see that the beat note frequency noise is approximately 
white (Fig.[^. 

The footprint and complexity of the lock is reduced significantly using this ap¬ 
proach as each additional slave laser merely requires a beam splitter to overlap the 
two laser beams, a focusing lens, and a photodiode, all of which can easily fit on a 
small breadboard with beams delivered via optical fibers. Additionally, the cost to 
implement the offset lock is more than an order of magnitude smaller than the cost 
of alternative approaches such as an additional transfer cavity or a hollow-cathode 
lamp. 
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Fig. 3: (a) Energy level diagram of a three level A system. Here, two of the hyperfine 
ground states of Yb^, separated by 12.6 GHz, are used as the qubit states |0) and 
11). (b) Optical frequency combs from a single mode-locked laser. The combs are 
shifted relative to each other to ensure that the ion is driven by a pair of comb teeth 
with a frequency difference equal to the qubit frequency /qubit- 


3 Digital optical frequency comb lock for Raman transitions 

For trapped ion qubits, driving qubit transitions using optical fields is convenient as 
they can be tightly focused for the individual addressing of single ions in a closely 
packed linear chain of ions. Optical frequency combs from mode-locked lasers are 
ideal for this purpose due to their broad bandwidth (10 GHz to 100 THz) which can 
bridge the wide energy separation of atomic states, and tightly spaced comb teeth, 
which are separated by the laser repetition rate (/rep), usually in the 50-1,000 MHz 
range nsni. 

The magnetically insensitive hyperfine ground states of the ^^^Yb^ ion serve as 
our qubit states with |0) =^ 81 / 2 !^ = 0,my = O) and |1) =^ 81 / 2 !^ = 1,^/ = O), at 
a qubit frequency /qubit of 12.6 GHz. This system can be approximated as a three 
level A system, shown in Fig. [^a). Two-photon stimulated Raman transitions be¬ 
tween qubit states | 0 ) and 11 ) can be driven using two optical frequency combs ifTSll . 
Fig. [^b) shows optical frequency combs generated from a single mode-locked laser 
where the relative offset between the two combs, /2 — /i, is introduced by an acousto¬ 
optic modulator (AOM). When the frequency difference /i — /2 is such that 

/qubit = nX /rep =b (/i - f2) (1) 

where n is an integer, an off-resonant two-photon Raman process induces the tran¬ 
sition between |0) and |1) states through the virtual excited state \e). Because the 
Raman transitions are driven by the comb teeth pairs separated by exactly /qubit, it is 
sufficient to stabilize the relevant pair of comb teeth to each other, instead of stabiliz¬ 
ing the absolute frequency or repetition rate of the laser output. 
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Thermal drifts and the vibration of mechanical parts in the laser cavity are among 
the main causes of drift in the laser repetition rate, which decreases gate fidelities 
by detuning the frequency comb pairs from /qubit- Stabilization of /rep can be accom¬ 
plished by directly stabilizing the laser cavity length; While this method provides suf¬ 
ficient stability and was used in previous experiments (TTl, the lock bandwidth is typ¬ 
ically limited to the kHz range due to its reliance on piezo-translation of a relatively 
large mirror mi. This method requires access to the laser cavity, not available in 
many commercial lasers used in these experiments, and requires a stable microwave 
oscillator that filters near /qubit which can be costly. Our alternative approach is to 
measure the drift of /rep by locking a digitally-controlled oscillator (such as a DDS) 
to the repetition rate signal using a digital phase-locked loop (PLL), and compute 
the driving frequency, / 2 , of the optical modulator to satisfy Eq. We implement 
the lock digitally, using a proportional-integral (PI) controller implemented in a field 
programmable gate array (FPGA). The system provides a low-cost and scalable al¬ 
ternative to hardware-intensive analog methods and is a viable alternative when laser 
cavity access is limited. 


3.1 Implementation 

The optical frequency comb stabilization setup is shown in Fig. A mode-locked 
titanium sapphire laser (centered at 752 nm, 2 ps FWHM pulse duration) is used to 
generate an optical frequency comb with comb teeth separation of /rep ~ 76 MHz. 
The laser output is frequency doubled via second harmonic generation (SHG) using 
a BiBO crystal (length 8 mm). The light from the SHG output (376 nm) is split into 
two approximately co-propagating frequency combs by driving an AOM with two 
frequencies: /i, controlled by the main experimental control FPGA and used to scan 
the detuning of the Raman beams in the experiments, and / 2 , used to stabilize the 
frequency difference between the pairs of comb teeth that drive the Raman transition. 
The co-propagating beam geometry ensures temporal overlap of the pulses at the ion. 

A small fraction of the laser output at the fundamental wavelength (752 nm) is 
incident on a photodiode (PD, S9055 Hamamatsu), and then amplified and low-pass 
filtered to attenuate the higher harmonics of /rep. The filtered PD output is mixed with 
the amplified output of a DDS (DDSO, Analog Devices, AD9912, 4 /iHz frequency 
resolution). The mixer output is low-pass filtered to eliminate the high frequency 
harmonics, while preserving the portion related to the frequency and phase difference, 
which provides the error signal for the PI controller in the PLL. Digitization of the 
error signal is performed by an analog-to-digital converter (ADC, Analog Devices 
AD7671, 16-bit resolution, 1 MSPS) and then read by the FPGA (Altera Cyclone II). 
The PI computation is performed in the FPGA which sets the frequency values of 
DDSO (/o) to track /rep, and DDS2 (/ 2 ) to satisfy Eq. All DDSs share the same 
reference clock. 

The lock operates as follows: a small change in /rep, in time step k, results in 
a digitized error signal ^ (/rep — /o)^^ + {^rep ~ where At the sampling 
interval, and (f^rep and 0o are the phases of the PD and DDSO outputs, respectively. 
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Fig. 4: Schematic of setup used to stabilize the frequency difference between the pairs 
of comb teeth that drive the Raman transition. Laser light is incident on a photodiode 
(PD), which is used, in conjunction with a mixer, to lock the frequency of DDSO (/o) 
to the repetition rate (/rep). The FPGA sets the frequency of DDS2 (/2) based on /o to 
keep the difference frequency of the relevant frequency combs, shifted by /i and / 2 , 
constant at /qubit- SHG: second harmonic generation, DDS: direct digital synthesizer, 
ADC: analog-to-digital converter, PD: photodiode, FPGA: field programmable gate 
array, AOM: acousto-optic modulator. 


Given constants P and /, each time step a correction is made to /o where 

k 

m+\) = Mk) + iPxek) + iIx'£ e„). (2) 

n=l 

This locks the frequency and phase of fo(t) to /rep(0 feedback control. Conse¬ 
quently, when /rep drifts such that 


/rep(^+ 1) ~ /rep(^) + ^5 (3) 

S is known, as the corresponding /o is set by the FPGA to match /rep. The FPGA 
feeds-forward the value of the harmonic of 5 to / 2 , so that 

f2(k+l)=f2(k)+nxS, (4) 

This maintains the equality in Eq. and 

/qubit = ^ X [/rep(^ -f 1) + 5] -f (/i — [f2{k) + ^ X 5]), (5) 

allowing for stable qubit transitions as /rep drifts. With /rep ~ 76 MHz and /qubit ~ 
12.6 GHz, n = 166. 
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Time (s) 

Fig. 5: Lock response to a finite /rep step. By replacing the PD in Fig. with a DDS, 
we can measure the reaction of the lock to a finite jump in /rep for the chosen P and I 
coefficients. We choose an I that, with P = 1, will produce a slew rate of ~ 50 Hz/s. 
The difference between the current and initial DDS frequency is shown as filled blue 
circles, and the open red circles represent the ADC value read. 


3.2 Lock parameters 

By examining the mixer output (error signal) on an oscilloscope, we determined a 
maximum rate of /rep drift of^5 Hz/s. To establish the optimal integral gain coeffi¬ 
cient, /, and the frequency difference detection sensitivity of our system, we replaced 
the PD output with another DDS, which was amplified to the same output power level 
and set to about the same frequency as the PD output. After enabling the lock with a 
proportional coefficient P = 1, the DDS frequency was changed abruptly, as shown in 
Fig. The coefficient I was chosen to provide a 50 Hz/s slew rate, sufficient for the 
lock to follow our expected 5 Hz/s rate of change. With the I determined, to minimize 
noise we used the minimum proportional coefficient P necessary to stabilize the lock. 


3.3 Digital filtering 

In our lock implementation, digital filters can be inserted without any hardware mod¬ 
ification. A simple averaging filter was implemented where, upon receiving N error 
signal samples at the maximum rate of 1 MSPS from the ADC, the PI loop proceeds 
using ek = as the error signal. Fig. j^a) shows the ADC readings where 

/rep is again replaced with a DDS. The DDS is set to take finite frequency steps, sim- 
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Fig. 6: Digital averaging. With the frequency comb stabilization system unlocked 
the PD is replaced with a DDS, set to take finite frequency steps, (top row) When 
/rep —fi 7^0 the error deviates from zero until /rep — /2 = 0 again, (bottom row) The 
error signal noise is reduced by digital averaging. Using a simple linear fit, we see 
values of 0.121, 0.034, 0.007 (V^) for A/^ =1, 4, and 16 respectively. 


ilar to Fig. with the system unlocked, at different values of N. The ADC noise is 
clearly reduced by digital averaging as evidenced by the residuals of a simple linear 
regression, shown in Fig.j^b). While averaging decreases the error signal noise, the 
locking bandwidth is also reduced. With no digital averaging, the locking bandwidth 
of the digital PI loop is limited to 1 MHz by the conversion speed of the ADC. With 
digital averaging, the lock bandwidth is reduced by 1/A. 


3.4 Lock verification 

High fidelity qubit operations require that Raman beat notes stay coherent with the 
qubit frequency for the duration of the computation. To measure the coherence time, 
a Ramsey experiment was performed. First, the qubit is initialized into the |0) state, 
followed by a ;r/2 pulse that rotates the qubit into a superposition state. After a wait 
time of T, a second Tljl pulse, with a phase shift of (j) with respect to the first pulse, 
is implemented followed by a measurement. This procedure is repeated for many 
values of (j) resulting in a fringe. The fringe visibility indicates how well difference 
frequency of the Raman beams remains synchronized with the phase of the qubit. 
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Ramsey Wait Time r (ms) 

Fig. 7: Coherence time using the digital optical frequency comb lock. Red data points 
represent the Ramsey fringe visibility at various values of time delay, T. The fringe 
visibility as a function of T is fit to a Gaussian function, black line, resulting in a 
coherence time of 811 ± 69 ms. Dark gray represents standard error on the fit. 


Using the digital optical frequency comb lock we measure a coherence time of 
811 ± 69 ms, shown in Fig.[^ The coherence time is determined by fitting the Ram¬ 
sey fringe visibility versus T to the Gaussian function Axexp(-T^/a^), where a is 
the coherence time. By directly stabilizing the cavity length in the same laser, we 
previously measured a coherence time of 687 ± 55 ms using a similar procedure (nl. 


3.5 Global synchronization 

In a scalable quantum hardware system, it is crucial that the phase coherence of all 
classical fields that manipulate various parts of the system be maintained. In the case 
of Raman beams that drive quantum logic gates, the phase of the Raman beat notes 
must be maintained, in case there are multiple pairs of Raman beams driving quan¬ 
tum gates across the system. We ensure this by driving all the DDS controllers and 
microwave sources in the system from a single clock signal derived from an atomic 
standard. The relative phase difference between all of the RF and microwave control 
signals reaching the qubits must be calibrated (and stabilized) to operate a distributed 
quantum system such as MUSIQC. 
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4 Digital optical intensity lock 

Intensity fluctuations of the Raman beams at the ion location result in the degradation 
of gate fldelities, and drifts in Rabi frequency and Stark shift calibrations. For high 
fldelity gates, control of the laser power at the ion location is critical. Here, we imple¬ 
ment a power stabilization system using a PD, ADC, FPGA, and a DDS. As shown 
in Fig.[^ after the frequency combs exit the flber, a small amount of collimated light 
is picked-off and placed on a PD. The PD signal is digitized by an ADC (AD7671) 
and relayed to the FPGA (Altera, Cyclone II). This provides the error signal for a 
PI feedback loop, similar to that used in Sect. |3.1[ The FPGA actively stabilizes the 
power on the PD by controlling the amplitude of the DDS signal driving AOM 1. This 
implementation allows for gating of the PI loop by a digital trigger pulse, sent from 
the FPGA in our main controller. In most experiments the Raman beams must be 
turned on and off. A digital trigger allows the PD to be placed after the AOMs, closer 
to the ion, providing better intensity stabilization. In our experiments, we turn the 
Raman beams and PI loop on during the Doppler cooling cycle, allowing the power 
in the Raman beams to be stabilized before the experiment. Once the cooling cycle is 
complete, the feedback control is turned off, and we hold the flnal value of the AOM 
amplitude setting from the lock throughout the remainder of the experiment until the 
next cooling cycle. Using this sample and hold approach, slow intensity fluctuations 
due to beam pointing instability before the flber are corrected. 

Figure]^ shows the impact of intensity stabilization. When the lock is turned off 
(constant DDS amplitude, as in Fig. |^), the power delivered through the flber drifts 
as monitored by the PD (sampled by the ADC). Once the locked is turned on, the 
DDS amplitude is adjusted to maintain constant optical power out of the flber (Fig. 

Wpy 


5 Next generation lock 


We are developing a general-purpose 8-channel digital proportional-integral-differential 
(PID) control system to replace the PI loop implementations discussed in sections 2-4 
and to be used for future locking tasks. It supports multiple user-selectable outputs 
for the lock, such as DC voltage through digital-to-analog converters (DACs) and 
frequency or amplitude of DDS signals, among others. 


The PID logic is implemented on a Xilinx Spartan 6 FPGA, detailed in Fig. 10 


The FPGA is mounted on a circuit board that also contains an 8-channel ADC (Ana¬ 
log Devices AD7608, 18-bit resolution, 200 ksps) and DACs (Texas Instruments 
DAC8568, 16-bit resolution) for each channel. The locking system can support eight 
PID locks concurrently, each with optional DC output, using the integrated DACs 
(66 kHz maximum update rate), or AC output, via an off-chip DDS (Analog Devices 
AD9912, 100 kHz maximum update rate). A Python graphical user interface (GUI) 
was designed for selection of locking parameters and real-time lock monitoring, with 
communication to the FPGA provided by a USB 2.0 link. Digital Altering, as in sec¬ 


tion |3^ is available on each channel through an oversampling module that performs 
digital averaging of the error signal. 
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Fig. 8: Experimental setup with optical intensity lock. A small fraction of the Raman 
laser light is incident on a photodiode (PD), which provides the error signal for a PI 
controller implemented in an FPGA. The FPGA dynamically sets the amplitude of 
the DDS signal driving AOM 1 to keep the intensity on the PD constant. 
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Fig. 9: Photodiode signal with and without digital optical intensity lock. The voltage 
read by the ADC (red crosses) is proportional to the optical intensity of the Raman 
beams, (a) When the intensity lock is off, the DDS amplitude (blue circles) is held 
constant and the ADC voltage, representing the optical power in the beam, drifts, 
(b) With the intensity lock on, the DDS amplitude changes to keep the ADC voltage 
constant. 
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Fig. 10: Next generation PID system. The ADC controller reads eight channels of er¬ 
ror data from the AD7608. The AD7608 runs at its maximum update rate of 200kHz. 
Error data is passed, one channel at a time, from the ADC controller into the PID 
processing pipeline. The first stage of the pipeline is an oversampling module. The 
module maintains a running average of data words as they are received. The aver¬ 
age is passed to the PID filter stage when the channel oversample ratio is satisfied. 
The PID filter applies a discrete PID filter to the averaged error data and the result 
is routed to the appropriate output channels. The output processor enforces signal 
bounds and optionally applies a linear transformation before passing the final data 
signal to the appropriate DAC or DDS controller. 


6 Digital to Analog Converter System 

Historically, macroscopic 4-rod RF-Paul ion traps have been used for trapped ion 
quantum information experiments Ill9ll20l[6ll . However, in recent years surface-electrode 
ion traps have come into use, as they provide repeatable trap parameters and support 
increased complexity. The trapping potential is constructed from a sum of RF and DC 
components. In a surface-electrode ion trap, these fields are generated by applying 
voltages to microfabricated metal electrodes located on a planar surface at a distance 
on the order of 100 jdm from the ions. The voltages applied to these electrodes can be 
modulated to shuttle the ions along the trap axis or split and re-combine ion chains 
with high precision EHEHa. While ion shuttling is central to many quantum com¬ 
puting schemes, this introduces the need for a multi-channel DC-voltage supply with 
precise timing control between channels. Although only a single RF source is re¬ 
quired for the RF trapping potential, typical surface traps require 40-100 DC voltages 
in order to accommodate the desired shuttling range and precise motional control of 
the ions. Any electric field noise near the motional frequencies of the trap (typically 
1 MHz to 5 MHz) will cause unwanted heating of the ion motion in the trap. Existing 
commercial systems are expensive, cumbersome to program, and generate noise near 
the trap frequencies due to their synchronous update rate of ~ 1 MS/s ll^ . 

In an effort to create a low-cost scalable voltage solution, we designed a DAC 
system with 25 DAC8734 chips, each with 4 channels (schematic shown in FigpT]). 
An onboard FPGA module, the Opal Kelly XEM6010-FX45, communicates serially 
with each of these 25 DAC chips and enables real-time control of the voltages. This 
system achieves a five-fold cost reduction compared to a similar National Instruments 
DAC system, while also giving much better control over the update rates. Figure [TT] 
also shows that this system can drive multiple 100-channel DAC systems, in case the 
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Fig. 11: A schematic overview of the 100-channel DAC system. The voltages are ma¬ 
nipulated in real time by an Opal Kelly XEM6010 FPGA module which is controlled 
by a python GUI via a universal serial bus (USB) connection. The FPGA output pins 
are directly connected to 25 DAC8734 chips using a serial peripheral interface (SPI)- 
like protocol where each chip outputs 4 analog voltages. This SPI communication 
limits the update rate of the DACs to 430 kHz. The analog voltages go through several 
stages of RC filters before they arrive at the trap electrodes. These filters, which have 
a cutoff frequency of 150-400 kHz, reduce the noise on the trap electrodes near the 
resonant frequencies of the trap (around 1-5 MHz). The Opal Kelly module includes 
128 MB of random access memory (RAM) which allows for multiple DAC output 
configurations with fast switching between configurations. An aluminum heat-sink is 
attached to the DACs to dissipate up to 2 W of power per chip. 


overall quantum system involves multiple surface trap chips distributed over multiple 
vacuum chambers. 

The DAC8734 chip was chosen primarily for its asynchronous nature. It has been 
observed that synchronous DACs whose update rates are near the ion trap frequencies 
cause ion heating ||2T]| due to peaks in the electric field noise spectrum at the update 
rate and its harmonics. Most commonly used two-qubit gate schemes in ion trapping 
systems use the ions’ shared motional state as a bus mode to generate entanglement 
between the ions. This requires that the motional quanta of the bus mode stays con¬ 
stant throughout duration of the gate, and so the heating rate will significantly impact 
the achievable gate fidelity. An asynchronous DAC has the advantage that updates 
can be paused while gates are being performed. In this way, the effect of the digi¬ 
tal noise arising from the DACs on the fidelity of the two-qubit gates in the system 
can be minimized. Alternatively, synchronous DACs that update at a rate much faster 
than the ion trap frequencies have been found to cause little heating 12^ . However, 
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Table 1: A table of the electrical characteristics of the DAC8934. 


DAC Resolution 

16-bits 

Max Output Voltage 

±10 V 

Max Update Rate 

430 kHz 

Response Time 

20 ns 

Digital Clock Frequency 

< 50 MHz 

Max Settling Time (Output changes from -10 V to 10 V) 

8 jUs 


these DACs are both expensive and require significant computational overhead, as 
FPGA clock frequencies are typically on the same order as the update rates. Asyn¬ 
chronous DACs can be updated at an arbitrarily low rate, well away from the ion trap 
frequencies. Other than ion shuttling, our DACs are typically not updated, effectively 
removing the noise associated with the DAC update rate during operations that are 
sensitive to motional heating. Table shows some of the electrical properties of the 
DAC8734. The 16-bit resolution enables control over the voltage down to ±300/i V, 
and the response time is a single clock cycle (20 ns). This response time is the time 
the DAC takes to respond to a latch command after an update has been sent (which 
takes 2000 ns, or 100 clock cycles). 

For many ion trapping applications, it is important to be able to control the out¬ 
put of these DACs in real-time. We make use of a commercial FPCA module (Opal 
Kelly XEM6010-LX45). Not only does this enable forward compatibility with fu¬ 
ture FPGAs, but it also simplifies the design process, as the FPGA module includes 
random access memory (RAM), which enables memory-intensive operations such as 
ion shuttling. To control the system as a whole, we make use of multiple layers of 
software and hardware (shown in Figure pT|). A Python program allows the user to 
generate a series of voltage sets to be executed. This program passes assembly code 
to the FPGA which is executed in real time and controls the serial outputs to the DAC 
chips. The voltage sets are sent from the PC to a combination of block memory on the 
FPGA (32 sets maximum) and RAM (8 M sets maximum) on the Opal Kelly module. 
When the system is not in its shuttling mode, the FPGA uploads one set of voltages 
to the DAC chips and then turns off all outbound clocks to reduce the noise on the 
DAC outputs. The FPGA interpolates between these sets at a user defined update rate 
(< 430 kHz) and step size. 


6.1 DAC noise analysis 

The outputs of the DACs pass through several layers of filtering before reaching the 
ion trap electrodes in an effort to reduce noise near the trap frequencies. These filters 
reduce both noise generated by the DAC system itself and noise picked up on the 
wires that connect the DAC system to the ion trap. It is common practice in surface 
electrode ion traps to use the electrodes powered by the DACs as an AC ground for 
the RF voltages used inside the trap. If this RF voltage has a weak ground path, it can 
radiate throughout the system and generate yet more noise. Within these constraints, 
we found it necessary to provide a good path to ground for the system at the filters 
immediately outside of the electrical vacuum feedthrough. The RC filters are housed 
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Fig. 12: A schematic of the filters used between the DAC system and an ion trap. 
Outside of the vacuum chamber, an in-line Pi-filter is connected to an RC filter on a 
custom PCB that is explicitly grounded. The last set of filters are on the ion trap chip 
in the vacuum chamber. 


in a Faraday cage machined out of a block of steel that is connected directly to both 
the trap ground and to a grounding strap in an effort to reduce the noise on the DAC 
outputs. The filters used have C = 470 pF and R = 1 ki2, with a cut off frequency 
/ ~ 340 kHz. This connector provides a low-resistance path to ground and acts to 
shield the lines after the filters from any additional noise. Commercial in-line filters 
(API Technologies 56-705-005 filters) with a cutoff frequency of 800 kHz were added 
to increase the filtering at high frequency and avoid any unwanted pickup near the 
self-resonant frequency of the capacitors used in the RC filters. Figure shows a 
circuit diagram of these filters. 

It is useful to evaluate the noise performance of the DAC system in two modes: 
synchronous mode, where the output voltages are updated at their maximum rate, 
and asynchronous mode, where the output voltages are updated once and left con¬ 
stant. In practice, the update rate can be anywhere between these two extremes, but 
synchronous mode is a good approximation of the behavior of the system during ion 
shuttling, and asynchronous mode is a good approximation for all other times. Figure 

shows the noise spectrum of the output of the DAC system for both of these modes. 
The noise spectrum for the synchronous case clearly shows noise at the update rate 
(430 kHz), its half-harmonic, and its harmonics. However, much of this noise will be 
significantly reduced by the filters shown in figure 

We have demonstrated the performance of a low noise DAC system, with lower 
cost compared to commercially available DACs. These DACs are optimized for the 
common case of static ions while retaining significant shuttling capabilities. Dur¬ 
ing shuttling their performance will be slower than the DACs shown by Bowler et. 
al. 1221 . yet comparable to the NI PXI-6733 used by Shu et. al. ||23l. Our DAC sys¬ 
tem is significantly less expensive, easier to program, and require less computational 
overhead than the other two alternatives. 
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Fig. 13: Noise measured on the a single output channel of the DAC system at 5 V. The 
spectrum analyzer resolution bandwidth was 100 Hz, and the input impedance was 
1 M^2. The system was run in two modes: ’’synchronous mode”, where the system 
is updated at the maximum rate of 430 kHz, and ’’asynchronous mode”, where the 
system is updated only once, (a) The synchronous mode noise spectrum is indicative 
of the system’s performance during the fastest possible shuttling operation, while (b) 
the asynchronous mode noise spectrum is analogous to the system’s behavior when 
not shuttling. These measurements do not include any filtering. 


7 Summary and conclusion 

As we scale up trapped ion quantum computing to a system large enough to perform 
significant quantum computations, the amount of classical hardware required will 
also grow. Here we described our efforts towards designing control hardware that is 
scalable and extensible while minimizing expense and spatial requirements. To this 
end, we have developed an offset frequency lock for CW lasers, such as those used for 
Doppler cooling, which allows many slave lasers to be locked to a single master laser, 
where each slave laser can also serve as a new master laser in a tree like structure. 
A single laser locked to a stable reference is sufficient to frequency stabilize a large 
number of slave lasers with little expense and space required for each slave laser lock. 
For pulsed laser used to drive Raman transitions, we have developed a pulsed laser 
Raman beam stabilization system that does not require access to the laser cavity, and 
can be implemented using simple and inexpensive low frequency components. With 
this lock the coherence time measured using a ^^^Yb^ ion is 811 (69) ms, which is 
sufficient a large scale computation containing thousands of logic gates. Currently we 
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are developing a general purpose PI loop system that could easily replace the pulsed 
laser Raman laser stabilization, intensity, and offset frequency locking controllers and 
software. The compact and economical system will be able to support 8 simultaneous 
locks with DC output, or optional off-chip DDS output. As each logic unit in the 
MUSIQC architecture will likely require an individual surface electrode ion trap with 
100 or more DC electrodes, we have developed a low-noise, asynchronous, DAC 
system that compares favorably to the available commercial systems in cost, size, 
and flexibility. 
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